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ABSTRACT 



Context. Understanding the 3D structure of coronal magnetic field is important to understand: the onset of flares and coronal mass 
ejections, the stability of active region, and to monitor the magnetic helicity & free magnetic energy and other phenomena in the solar 
atmosphere. Routine measurements of the solar magnetic field are mainly carried out in the photosphere. Therefore, one has to infer 
the field strength in the higher layers of the solar atmosphere from the measured photospheric field based on the assumption that the 
corona is force-free. Meanwhile, those measured data are inconsistent with the above force-free assumption. Therefore, one has to 
apply some transformations to these data before nonlinear force-free extrapolation codes can be applied. 

Aims. Extrapolation codes in cartesian geometry for modelling the magnetic field in the corona do not take the curvature of the Sun's 
surface into account and can only be applied to relatively small areas, e.g., a single active region. Here we apply a method for nonlinear 
force-free coronal magnetic field modelling and preprocessing of photospheric vector magnetograms in spherical geometry using the 
optimization procedure. 

Methods. We solve the nonlinear force-free field equations by minimizing a functional in spherical coordinates over a restricted area 
of the Sun. We extend the functional by an additional term, which allows to incorporate measurement error and treat regions with 
lacking observational data. We use vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey 
(SOLIS) to model the coronal magnetic field. We study two neighbouring magnetically connected active regions observed on May 15 
2009. 

Results. For vector magnetograms with variable measurement precision and randomly scattered data gaps (e.g., SOLIS/VSM) the 
new code yields field models which satisfy the solenoidal and force-free condition significantly better as it allows deviations between 
the extrapolated boundary field and observed boundary data within measurement errors. Data gaps are assigned to an infinite error. 
We extend this new scheme to spherical geometry and apply it for the first time to real data. 

Key words. Magnetic fields - Sun: corona - Sun: photosphere - methods: numerical 
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1. Introduction 

Observations have shown that physical conditions in the solar 
atmosphere are strongly controlled by solar magnetic field. The 
magnetic field also provides the link between different manifes- 
tations of solar activity like, for instance, sunspots, filaments, 
flares, or coronal mass ejections. Therefore, the information 
about the 3D structure of magnetic field vector throughout the 
solar atmosphere is crucially important. Routine measurements 
of the solar vector magnetic field are mainly carried out in the 
photosphere. Therefore, one has use numerical modelling to in- 
fer the field strength in the higher layers of the solar atmosphere 
from the measured photospheric field based on the assumption 
that the corona is force-free. Due to the low valu e of the plas ma /? 
(the ratio of gas pressure to magnetic pressure) dGarvll200ll) . the 
solar corona is magnetically dominated. To describe the equilib- 
rium structure of the coronal magnetic field when non-magnetic 
forces are negligible, the force-free assumption is then appropri- 
ate: 



(VxB)xB = 



V-B = 
B = H b s on photosphere 



(2) 
(3) 



where B is the magnetic field and H b S is 2D observed surface 
magnetic field on photosphere. Extrapolation methods have been 
developed for different types of force-free fi elds: potential field 
extrapolation dSchmidtj fl96l ISemell fl967h. linear force-free 
field extrapol ation dChiu & Hiltonl 1 1 977b ISeehafedll978l Il982b 
ISemell 11988b IClegg et al.l l2000h. and nonlinear force-free field 
extrap o lation dSakurail 1981b IWuetaL 1990b Cuperman et al.l 



1991 



Demoulin et al 



lRoumeliotisril996l 



1992F 



Mikic & McClvmont 1994; 
lAmarietal.l 1 19971 1 1999b lYan & Sakurail 
| 2000b IValori et al ll2005b IWheatlandl 12004b IWiegelmannl l2004b 
Amari et alj2006bllnhester & Wiegelmann|2006l) . Among these, 
the nonlinear force-free field has the most realistic description 
of the coronal magnetic field. For a more complete review of 
existing methods for computing nonlinear force-free coronal 



magne ti c fields, we refe r to the r eview papers b y Amari et all 
>97l). ISchriiver et all d2006l) . iMetcalf et alj d2008l) . 



(1997) 



(1) IWiegelmannl(l2008h 



and 
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The magnetic field is not force-free in the photosphere, 
but becomes force-f ree roughly 400 km above the photosphere 
dMetcalf et al.l 1 19951) . Nonlinear force-free extrapolation codes 
can be applied only to low plasma-/? regions, where the force- 
free assumption is justified. The preprocessing scheme as used 
until now modifies observed photospheric vector magnetograms 
with the aim of approximating the magnetic field vector at 
the bottom of the force-free domain (|Wiegelmann etldl 120061 : 
iFuhrmann et alj2007t iTadesse et al.l2009l) . The resulting bound- 
ary values are expected to be more suitable for an extrapolation 
into a force-free field than the original values. Preprocessing is 
important for those NLFF-codes which uses the magnetic field 
vector on the boundary directly. Consistent computations for the 
Grad-Rubin method, which uses B„ an d J„ (or a) as bound- 
ary co ndition have been carried out by IWheatland & Regnierl 
(120091) . 

In this paper, we use a larger computational domain which 
accommodates most of the connectivity within the coronal re- 
gion. We also take uncertainties of measur ements in vector mag - 
netograms into accounts as suggested in iDeRosa et al.l (120091). 
We im plement a preprocessing procedure of ITadesse et al.l 
(120091) to SOLIS data in spherical geometry by considering the 
curvature of the Sun's surface for the large field of view contain- 
ing two active regions. We use a spherical version of the opti- 
mization procedure that has been implem ented in cartesian ge- 
ometry in lWiegelmann & Inhesterld2010h for synthetic boundary 
data. 



2. Method 

2.1. The SOLIS/VSM instrument 

In this study, we use vector magnetogram observations from 
the Vector Spectromagnetograph (VSM; see iJones et al.ll2002l) . 
which is part of the Synoptic Optical Long-term I nvestigations 
of the Sun (SOLIS) synoptic facility (SOLIS; see iKeller et al.l 
I2003I) . VSM/SOLIS currently operates at the Kitt Peak National 
Observatory, Arizona, and it has provided magnetic field obser- 
vations of the Sun almost continuously since August 2003. 

VSM is a full disk Stokes Polarimeter. As part of daily syn- 
optic observations, it takes four different observations in three 
spectral lines: Stokes /(intensity), V (circular polarization, Q, 
and U (linear polarization) in photospheric spectral lines Fe I 
630.15 nm and Fe I 630.25 nm , Stokes / and V in Fe I 630.15 
nm and Fe I 630.25 nm, similar observations in chromospheric 
spectral line Ca II 854.2 nm, and Stokes / in the He I 1083.0 
nm line and the near-by Si I spectral line. Observations of /, Q, 
U, and V are used to construct a full disk vector magnetograms, 
while I-V observations are employed to create separate full disk 
longitudinal magnetograms in the photosphere and the chromo- 
sphere. 

In this study, we use a vector magnetogram observed on 15 
May 2009. The data were taken with 1.125 arcsec pixel size and 
2.71 pm spectral sampling. ( In December 2009, SOLIS/VSM 
has upgraded its cameras from Rockwell ( 90 Hz, 18 micron 
pixels ) to Sarnoff ( 300 Hz, 16 micron pixels ). This camera 
upgrade has resulted in improved spatial and spectral sampling 
). The noise level for line-of-sight component is about 1 Gauss. 
However, noise due to atmospheric seeing may be much larger, 
and the final measurement error depends on the measured flux, 
its spatial distribution as well as the seeing conditions. A rough 
estimate suggests a noise level of a few tens of Gauss for areas 
with a strong horizontal gradient of magnetic field and about 1 
arcsec atmospheric seeing. 



To create a single magnetogram, solar disk image is scanned 
from terrestrial South to North; it takes about 20 minutes to 
complete one vector magnetogram. After the scan is done, the 
data are sent to an automatic data reduction pipeline that in- 
cludes dark and flat field correction. Once the spectra are prop- 
erly calibrated, full disk vector (magnetic field strength, incli- 
nation, and azimuth) magnetograms are created using two dif- 
ferent approaches. Quick-loo k (QL) vec t or ma gnetograms are 
created based on algorithm bv lAuer et al.l (119771) . The algorithm 
uses the Milne-Eddington model of solar atmosphere, which 
assumes that magnetic field is unifor m (no gradie nts) through 
the layer of spectral line formation (lUnnol 1956b . It also as- 
sumes symmetric line profiles, disregards magneto-optical ef- 
fects (e.g., Faraday rotation), and does not separate contribution 
of magnetic and non-magnetic components in spectral line pro- 
file (i.e., magnetic filling factor is set to unity). A complete in- 
ve rsion of spectral data is don e later using technique developed 
bv ISkumanich & Litesl (119871) . This latter inversion (called ME 
magnetogram) also employs Milne-Eddington model of atmo- 
sphere, but it solves for magneto-optical effects and determines 
magnetic filling factor (fractional contribution of magnetic and 
non-magnetic components to each pixel). The ME inversion is 
only performed for pixels with spectral line profiles above the 
noise level. For pixels below the polarimetric noise threshold, 
magnetic field parameters are set to zero. 

From the measurements the azimuths of transverse magnetic 
field can only be determined with 180-degree ambiguity. This 
ambiguity i s resolved using th e Non-Potential Field Calculation 
(NPFC; see lGeorgoulisll2005l) . The NPFC method was selected 
on the basis of comparative investi gation of several m ethods 
for 180-degree ambiguity resolution dMetcalf et al.ll2006l) . Both 
QL and ME magnetograms can be used for potential and/or 
force-free field extrapolation. However, in strong fields inside 
sunspots, the QL field strengths may exhibit erroneous decrease 
inside sunspot umbra due to, so called magnetic saturation. For 
this study we choose to use fully inverted ME magnetograms. 
Fig. Q] shows a map of the radial component of the field as a 
contour plot with the transverse magnetic field depicted as black 
arrows. For this particular dataset, about 80% of the data pixels 
are undetermined and as a result the ratio of data gaps to total 
number of pixels is large. 

2.2. Preprocessing of SOLIS data 

The preprocessing scheme o flTadesse et ail (120091) involves min- 
imizing a two-dimensional functional of quadratic form in spher- 
ical geometry as following: 



H = argtnin(L p ) 

L p = HlL\ + H2L2 + yU 3 L3 + 



(4) 



where H is preprocessed surface magnetic field from the input 
observed field H oos . Each of the constraints L n is weighted by a 
yet undetermined factor fi„. The first term (n = 1) corresponds to 
the force-balance condition, the next (n = 2) to the torque-free 
condition, and the last term (n = 4) controls the smoothi ng. The 
explic it form of L\, L2, and L4 can be found in ITadesse et al.l 
(120091) . The term (n = 3) ensures that the optimized boundary 
condition agrees with the measured photospheric data. In the 
case o f SOLIS/VSM data we modified L3 with respect to the 
one in Tadess e et al 1 (120091) as follows, to treat those data gaps. 



W(0,<f>)-(H-H obs ), 



(5) 
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Fig. 1. Surface contour plot of radial magnetic field vector and vector field plot of transverse field with black arrows. 



In this integral, W(9, 0) = diag(w m ^ a i, Wtrans, w aans ) is a diagonal 
matrix which gives different weights to the different observed 
surface field components depending on their relative measure- 
ment accuracy. A careful choice of the preprocessing parameters 
fi„ ensures that the preprocessed magnetic field H does not devi- 
ate from the original observed field H b S by more than the mea- 
surement errors. As the result of parameter study in this work, 
we found fi\ = fi2 = 1-0, fi3 = 0.03 and //4 = 0.45 as optimal 
value. 

2.3. Optimization principle 

Equations (Q]i and (f2| can be solv ed with the help of an opti- 
mization princ iple, as proposed by IWheatland et al.l (l2000t) and 
generalized by IWiegelma nn ( 2004) for cartesian geometry. The 
method minimizes a joint measure of the normalized Lorentz 
forces and the divergence of the field throughout the volume 
of interest, V. Throughout this minimization, the photospheric 
boundary of the model field B is exactly matched to the observed 
H„bs and possibly preprocessed magnetogram values H. Here, 
we use the optimization approach for functional (L w ) in spherical 
geometry dWiegelmannl 120071: iTadesse et alj 120091) along with 
the new method which instead of an exact match enforces a min- 
imal deviations between the photospheric boundary of the model 
field B and the magnetogra m field H ODS by adding an appro priate 
surface integral term L p i loto (IWiegelmann & Inhestertl2010l) . 

B = argmin(L w ) 



La> - Lf + Ld + vLpko t0 (6) 

Lf = I 0) f (r, 9, (f>)B~ 2 \(V xB)xB\ 2 r 2 sin 0drd8d<p 
Jv 

L d = I u) d {r, 9, 0)| V • B\ 2 r 2 sin 8drd9d(f> 
Jv 

Lphoto = f (B - H ohs ) ■ W(6, <P)-{B- H obs )r 2 sin 6d6dcf> 

where Lf and L ( / measure how well the force-free Eqs. (Q} 
and divergence-free (f2]i conditions are fulfilled, respectively. 
ojf(r, 9, <p) and <wj(r, 9, <p) are weighting functions. The third inte- 
gral, Lpi lo t , is surface integral over the photosphere which allows 
us to relax the field on the photosphere towards force-free solu- 
tion without to much deviation from the original surface field 
data. W(9, 0) is the diagonal matrix in Eq. (0. 

Nume rical tests of the effect of the new term L p i wto were per- 
formed bv lWiegelmann & Inhestej(l2010h in cartesian geometry 
for syn thetic magnetic fie ld vector generated from Low & Lou 
model (lLow & Loul 19901) . They showed that this new method to 
incorporate the observed boundary field allows to cope with data 
gaps as they are present in SOLIS and other vector magnetogram 
data. Within this work, we use a spherical geometry for the full 
disk data from SOLIS. We use a spherical grid r, 9, <p with n r , rig, 
tiif, grid points in the direction of radius, latitude, and longitude, 
respectively. The method works as follows: 

- We compute an initial source surface potential field in the 
computational domain from H ro i, s , the normal component of 
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the surface field at the photosphere at r — IRq- The compu- 
tation is performed by assuming that a currentless (7 = 
or V x B = 0) approximation holds between the photo- 
sphere and some spherical surface S s (source surface where 
the magnetic field vector is assumed radial). We computed 
the solution of this boundary-value problem in a standard 
form of harmonic expansion in terms of eigen-solutions of 
the Laplace equation written in a spherical coordinate sys- 
tem, (r, 9, <f>). 

- We minimize L w (Eqs. [6]i iteratively without constraining 
H b s at the photospher e boundary as in prev ious version 
of Wheatland algorithm (I Wheatland et al.ll2000t) . The model 
magnetic field B at the surface is gradually driven towards 
the observations while the field in the volume V relaxes to 
force-free. If the observed field is inconsistent, the difference 
B - H„b s or B - H (for preprocessed data) remains finite de- 
pending in the control parameter v. At data gaps in H b S , we 
set Wradiai = and w, nm s = and respective field value is 
automatically ignored. 

- The state L w — corresponds to a perfect force-free and 
divergence-free state and exact agreement of the boundary 
values B with observations H i,s in regions where w raC iiai and 
Wtrans are greater than zero. For inconsistent boundary data 
the force-free and solenoidal conditions can still be fulfilled, 
but the surface term L p i wto will remain finite. This results in 
some deviation of the bottom boundary data from the ob- 
servations, especially in regions where w ra diai an d w tra „s are 
small. The parameter v is tuned so that these deviations do 
not exceed the local estimated measurement error. 

- The iteration stops when L w becomes stationary as 
AL W /L tt < 1(T 4 . 



3. Results 

We use the vector magnetograph data from the Synoptic Optical 
Long-term Investigations of the Sun survey (SOLIS) to model 
the coronal magnetic field. We extrapolate by means of Eq. (O 
both the observed field H a bs measured above two active re- 
gions observed on May 15 2009 and preprocessed surface field 
(H, that obtained from H i,s applying our preprocessing proce- 
dure). We compute 3D magnetic field in a wedge-shaped com- 
putational box V, which includes an inner physical domain V 
and the buffer zone (the region outside the physical domain), 
as shown in Fig. [3] of the bottom boundary on the photosphere. 
The wedge-shaped physical domain V has its latitudinal bound- 

o o 

aries at # m ; n = 3 and 9 mRX = 42 , longitudinal boundaries at 

o o 

<pmin = 153 and <p max = 212 , and radial boundaries at the pho- 
tosphere (r = IRq) and r = 1.75/?©. 

The weighting function cjf and cjj in Ly and Ld in Eq. (0 
are chosen to be unity within the inner physical domain V 
and decline with a cos i ne profile in the b uffer boundary re- 
gion dWiegelmannl2004l ; lTadesse et alj|2009l) . They reach a zero 
value at the boundary of the outer volume V. The distance be- 
tween the boundaries of V and V is chosen to be nd = 10 grid 
points wide. The framed region in Figs.[3](a-i) corresponds to the 
lower boundary of the physical domain V with a resolution of 
132 x 196 pixels in the photosphere. The original full disc vector 
magnetogram has a resolution of 1788 x 1788 pixels out of which 
we extracted 142 x 206 pixels for the lower boundary of the com- 
putational domain V, which corresponds to 550Mra x 12QMm on 
the photosphere. 

The main reason for the implementation of the new term 
Lphoto in Eq. © is that we need to deal with boundary 



data of different noise levels and qualities or even lack some 
data points completely. SOLIS/VSM provides full-disk vector- 
magnetograms, but for some individual pixels the inversion from 
line profiles to field values may not have been successful inverted 
and field data there will be missing for these pixels. Since the old 
code without the L p i lo to term requires complete boundary infor- 
mation, it can not be applied to this set of SOLIS/VSM data. In 
our new code, these data gaps are treated by setting W = for 
these pixels in Eqs. ((6]). For those pixels, for which H b s was 
successfully inverted, we allow deviations between the model 
field B and the input fields either observed H i,s or preprocessed 
surface field H using Eqs. (|6]l and so that the model field can be 
iterated closer to a force-free solution even if the observations 
are inconsistent. This balan ce is controlled by the Lagran gian 
multiplier v as explained in IWiegelmann & Inhesterl d2010h . In 
this work we used w ra diai - 100if mmj for the surface fields both 
from data with preprocessing and without. 

Figure|2] shows the position of the active region on the solar 
disk both for SOLIS full-disk magnetogram Q SOHO/EIT im- 
age of the Sun observed at 195 A on the same day at 16:OOUT0 
As stated in section 2.3, the potential field is used as initial con- 
dition for iterative minimization required in Eq. [6] The respec- 
tive potential field is shown in the rightmost panel of in Fig. [2] 
During the iteration, the code forces the photospheric boundary 
of B towards observed field values H bs or H (for preprocessed 
data) and ignores data gaps in the magnetogram. A deviation 
between surface vector field from model B and H i,s or H (for 
preprocessed data) occurs where H b S is not consistent with a 
force-free field. In this sense, the term L p / !0(0 in Eq. © acts on 
Hgbs similarly as the preprocessing, it generates a surface field 
B instead of H from H b s which is close to H b s , but consistent 
with a force-free field above the surface. In Fig. [3] we therefore 
compare the option of the preprocessing and the new extrapo- 
lation code (Eq. [6J on H i,s- The figure shows the surface mag- 
netic field differences of the preprocessed, un-preprocessed and 
the potential surface fields. 

In order to deternime the similarity of vector components on 
the bottom surface, we calculate their pixel-wise correlations. 
The correlation were calculated from: 



vec 1/2 

where i>, and m, are the vectors at each grid point ; on the bottom 
surface. If the vector fields are identical, then C vec = 1; if D; ± u -, 
, then C vec = 0. Table [TJ shows correlations of the surface fields 
from B pre - Hobs (where B pre is the model field obtained from 
preprocessed surface field H using Eq. ©) and B mpre - H bs 
(where B unpre is the model field obtained from observed surface 
field Hobs using Eq. (O). We have computed the vector correla- 
tions of the two surface vector fields for the three components 
at each grid points to compare how well they are aligned along 
each directions. From those values in Table [TJ one can see that 
the preprocessing and extrapolation with Eq. © act on H a b s in 
a similar way. In Figj4] we plot magnetic field lines for the three 
configurations in which the vector correlations of potential field 
lines in 3D box to both the extrapolated NLFF with and without 
preprocessing data are 0.741 and 0.793, respectively. 

To understand the physics of solar flares, including the lo- 
cal reorganization of the magnetic field and the acceleration of 
energetic particles, one has to estimate the free magnetic energy 

1 http://solis.nso.edu/solis data.html 

2 http://sohowww.nascom.nasa.gov/data/archive 
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Fig. 2. Left: Full disc vector magnetogram of May 15 2009 at 16:02UT. Middle: SOHO/EIT image of the Sun on the same day 
at 16:00UT. Right: potential magnetic field line plot of SOLIS vector magnetogram atl6:02UT, that has been computed from the 
observed radial component. 
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Fig. 3. Top row: Radial surface vector field difference of a), modelled B without preprocessing and H bs b). modelled B pre and 
H b s c). initial potential and H„i, s . Middle row:Latitudinal surface vector field difference of d). modelled B without preprocessing 
and H bs e )- modelled B pre and H b S f). initial potential and H bs- Bottom row:Longitudinal surface vector field difference 
of g). modelled B without preprocessing and H b S h). modelled B pre and H bs 0- initial potential and H bs- The vertical and 
horizontal axes show latitude, 6 and longitude, <p in degree on the photosphere respectively. 
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(a) Potential field 



(b) Field from data before preprocessing 



(c) Field from preprocessed data 



Fig. 4. a) Some field lines for the Potential field reconstruction, b) Nonlinear force-free reconstruction from SOLIS data without 
preprocessing, c) Nonlinear force-free reconstruction from preprocessed SOLIS data. The panels show the same FOV as in Fig. [2] 
(right panel). 



Table 1. The correlations between the components of surface 
fields from (B pre - H ohs ) and (B""'' re - H obs ). 



V 






u 


C V ec 


^Qunpre 


H bs)r 


(B'"' e 


- H bs)r 


0.930 




H bs) 9 


(B pre 


— H B bs)e 


0.897 


(g'lnpre _ 


H B bs)(f, 


(B" re 


- Hobs)<f, 


0.875 



Table 2. The magnetic energy associated with extrapolated 
NLFF field configurations with and without preprocessing. 



Table 3. The currents and average a calculated from those pix- 
els which are magnetically connected. The currents are given in 
Ampere (A). 





Inside left 


Between left 


Inside right 




active region 


and right ARs 


active region 


Magnetic flux (I0 lv Gem 2 ) 


3.32 


4.61 


2.08 


Total current (10 s A) 


49.6 


1.58 


32.17 


Average a (MmT 1 ) 


2.49 


0.08 


1.62 



Model 



E„ lff (10 32 erg) £ free (10 32 e rg) 



No preprocessing 
Preprocessed 



37.456 
37.341 



4.915 
4.800 



available for such phenomena, this is the free energy which can 
be converted into kinetic and thermal energy. From the energy 
budget and the observed ma gnetic activity in the ac tive region, 
iRegnier & Priest! d2007al) and lThalmann et al.l d2008h have inves- 
tigated the free energy above the minimum-energy state for the 
flare process. We estimate the free magnetic energy from the 
difference of the extrapolated force-free fields and the potential 
field with the same normal boundary conditions in the photo- 
sphere. We therefore estimate the upper limit to the free mag- 
netic energy associated with coronal currents of the form 

£free = ^ j (B 2 nlff - B^r 1 sinOdrdOd^ (8) 



where B pot and B,,/// represent the potential and NLFF magnetic 
field, respectively. The free energy is about 5 x 10 32 e rg. The 
magnetic energy associated with the potential field configura- 
tion is found to be 32.541 x 10 32 e rg. Hence £«;// exceeds E po , 
by only 15%. Table|2]shows the magnetic energy associated with 
extrapolated NLFF field configurations with and without prepro- 
cessing. The magnetic energy of the NLFF field configuration 
obtained from the data without preprocessing is a slightly large 
from the preprocessed boundary field, as the preprocessing pro- 
cedure removes small scale structures. 

The electric current density calculated from Ampere's law, 
J — V x B/4n, on the basis of spatially sampled transverse 
magnetic fields varies widely over an active region. In order 



to investigate how errors in the vector magnetograph measure- 
ments produce er rors in the vertical electric current densities, 
iLiang et all (120091) have numerically simulated the effects of ran- 
dom noise on a standard photospheric magnetic configuration 
produced by electric currents satisfying the force-free field con- 
ditions. Even if the current density can be estimated on the pho- 
tosphere, it is not intuitively clear how the change in the current 
density distribution affe cts a coronal magnetic configuration. 
iRegnier & Priestl (l2007bl) studied such modifications in terms of 
the geometry of field lines, the storage of magnetic energy and 
the amount of magnetic helicity. FigO shows Iso-surface plots 
of current density above the volume of the active region. There 
are strong current configurations above each active regions. This 
becomes clear if we compare the total current in between each 
active region with the current from the left to the right active re- 
gion. These currents were added up from the surface normal cur- 
rents emanating from those pixels which are magnetically con- 
nected inside or across active regions respectively. The result is 
shown in Table [3] The active regions share a decent amount of 
magnetic flux compared to their internal flux from one polarity 
to the other. In terms of the electric current they are much more 
isolated. The ratio of shared to the intrinsic magnetic flux is or- 
der of unity, while for the electric current those ratios are much 
less, 1.58/49.6 and 1.58/32.17, respectively. Similarly we can 
calculate the average value of a on the field lines with the re- 
spective magnetic connectivity. The averages are shown in the 
second row of Table [3] The two active regions are magnetically 
connected but much less by electric currents. 
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Fig. 5. Iso-surfaces (ISs) of the absolute current density vector 
|/| = 100mA ■ nC 1 computed above the active regions. 

4. Conclusion and outlook 

We have investigated the coronal magnetic field associated with 
the AR 11017 on 2009 May 15 and neighbouring active region 
by analysing SOLIS/VSM data. We have used the optimization 
method for the reconstruction of nonlinear force-free coronal 
magnetic fields in sphe rical geometry by restricting the code t o 
limited parts of the Sun (i Wiegelmannl2007llTadesse et al.l2009h . 
Different from previous implementations our new code allows us 
to deal with lacking data and regions with poor signal-to-noise 
ratio in the extrapolation in a systematic manner because it pro- 
duces a field which is closer to a force-free and divergence-free 
field and tries to match the boundary only wher e it has been re- 
liably measured (IWiegelmann & Inhesterll20loh . 

For vector magnetograms with lacking data point where zero 
values have been replaced for the signal below certain threshold 
value, the new code relaxes the boundary and allows to fulfill 
the solenoidal and force-free condition significantly better as it 
allows deviations between the extrapolated boundary field and 
inconsistent observed boundary data. With the new term L p i ml0 
extrapolation from H i, s and H yield slightly the same 3D field. 
However, in the latter case the iteration to minimize Eq. © sat- 
urates in fewer iteration steps. At the same time, preprocessing 
does not affect the overall configuration of magnetic field and its 
total energy content. 

We plan to use this newly developed code for upcoming data 
from SDO (Solar Dynamics Observatory )/HMI ( Helioseismic 
and Magnetic Imager) when full disc magnetogram data become 
available. 
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